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Abstract 

The  analytic  solution  to  the  shallow  water  equations  in  a  closed  equatorial  basin  subject 
to  arbitrary  wind  forcing  is  presented.  Two  particular  examples,  including  the  transient 
response  of  the  basin  to  a  moving  “top-hat”  wind  pulse  representative  of  the  atmospheric 
Southern  Oscillation,  are  worked  out  in  detail.  The  analytic  results  compare  favorably 
with  previously  published  numerical  model  results  and  offer  new  potential  for  interpreting 
equatorial  ocean  variability-both  modeled  and  observed. 


List  of  Figures 


1  Zonal  current  anomaly,  a)  t  =  5,  contour  interval  =  .002  m  s-1;  b)  t  =  20, 

contour  interval  =  .02  m  s-1.  (Westward  currents  are  dashed.) .  13 

2  Zonal  current  anomaly,  a)  t  —  50,  contour  interval  =  .lms"1;b)(  =  100, 

contour  interval  =  .2  m  s-1 .  14 

3  uo(*,x)  vs  x  (Vertical  scale  is  m  s-1x  100) . 18 

4  u2(t,x)  vs  x  (Vertical  scale  is  m  s-1x  100) . 19 

5  u4(t,  x)  vs  x  (Vertical  scale  is  m  s-1x  100) . 20 

6  u6(t,  x)  vs  x  (Vertical  scale  is  m  s-1x  100) . 21 

7  tig (t,x)  vs  x  (Vertical  scale  is  m  s-1x  100) . 22 

8  Zonal  current  anomaly,  a)  t  =  5,  contour  interval  =  .001  m  s-1;  b)(  =  50, 

contour  interval  =  .2  m  s'1 .  24 

9  Zonal  current  anomaly,  a)  t  —  100,  contour  interval  =  .1  m  s_1;  b)t  =  200, 

contour  interval  =  .5  m  s-1 .  25 

10  Zonal  current  anomaly,  a)  t  =  400,  contour  interval  =  .1  m  s-1;  b)f  =  500, 

contour  interval  =  .5  m  s_1 .  27 


List  of  Tables 


1  Values  of  constants 


Acknowledgements 


This  work  was  supported  by  the  National  Science  Foundation  under  grant  ATM85- 
14012.  Computations  were  made  at  the  W.  R.  Church  computer  center.  Computer  graphics 
were  produced  by  Mrs.  O.  Haney. 


1 


1  Introduction 


The  collective  atmosphere-ocean  anomaly  known  as  the  El  Nino-Southern  Oscillation 
(ENSO)  has  received  a  great  deal  of  attention  in  recent  years.  One  common  framework 
within  which  the  system  has  been  studied  is  the  shallow  water  equations.  With  this  ap¬ 
proach,  the  atmosphere  and/or  the  ocean  is  treated  as  the  first  baroclinic  mode  of  the 
shallow  water  equations.  The  equations  governing  the  two  systems  are  the  same,  except 
for  their  equivalent  depths  (or  phase  speeds)  and  their  forcing  functions.  The  atmosphere 
is  generally  forced  by  a  source  term  in  the  continuity  equation,  representing  heating  at 
middle  to  upper  levels  of  the  atmosphere;  the  ocean  is  usually  forced  by  a  momentum 
source,  representing  surface  wind  stress. 

Shallow  water  models  of  an  equatorial  ocean  basin  have  been  solved  numeri¬ 
cally  by  Hurlburt  et.  al.(l976),  Busalacchi  and  O’Brien  (1981),  Busalacchi,  et.  al.  (1983), 
Weisburg  and  T^ug  (1983),  and  others,  using  both  hypothetical  and  observed  forc¬ 
ing  functions  to  represent  enhanced  surface  wind  stress  due  to  anomalous  westerly 
wind  events.  McCreary  (1976),  Cane  and  Sarachik  (1976),  Cane  and  Sarachik  (1981), 
Cane  (1979),  McCreary  and  Lukas  (1986),  and  others  obtained  analytic  solutions  to  sim¬ 
ilar  sets  of  equations.  Similar  techniques  have  been  used  by  Gill  (1980),  among  others, 
to  investigate  the  response  of  the  first  baroclinic  mode  of  the  global  atmosphere  to  heat¬ 
ing  due  to  enhanced  convection  associated  with  sea  surface  temperature  anomalies  in  the 
equatorial  Pacific. 
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An  obvious  extension  of  these  studies  has  involved  the  coupling  of  the  atmospheric  and 
oceanic  shallow  water  models  through  the  introduction  of  interactive  forcing  functions.  The 
resulting  interactive  models  have  more  degrees  of  freedom  than  the  individual  component 
models,  and  they  provide  additional  insight  regarding  possible  interdependencies  and  res¬ 
onances  of  the  coupled  ocean-atmosphere  system.  Normal  mode  solutions  to  such  coupled 
(but  laterally  unbounded)  systems  were  studied  by  Lau  (1981),  Philander  et.  al.(1984), 
Yamagata  (1985),  Rennick  and  Haney  (1986),  and  Hirst  (1986),  all  of  whom  found  the 
possibility  of  unstable  interactions  under  different  sets  of  somewhat  restrictive  assumptions. 
Numerical  solutions  to  the  coupled  equations  have  been  studied  under  fewer  restrictions 
by  Rennick  (1983),  Schopf  and  Suarez  (1988),  and  Hirst  (1988). 

One  feature  of  most  of  the  analytic  studies  cited  above  is  that  they  failed  to  account 
for  the  effects  of  reflections  at  the  eastern  and  western  boundaries  of  the  ocean  ba«in. 
The  ocean  response  study  of  McCreary  and  Lukas  (1986),  and  all  the  normal  mode  stud¬ 
ies  of  the  coupled  system  noted  above,  effectively  treated  an  unbounded  equatorial  ocean 
band,  with  no  continents.  However,  the  numerical  studies,  as  well  as  thecretical  work  by 
McCalpin  (1988),  have  shown  that  boundary  reflections  play  an  impri,'*;i'’t  role  in  par¬ 
titioning  energy  among  Kelvin  and  Rossby  modes,  and  thereby  affecting  the  structure 
of  disturbances,  and  the  energy  available  to  the  unstable  modes.  In  addition,  observa¬ 
tional  evidence  presented  by  White,  et.  al.  (1985)  suggests  that  reflection  from  the  western 
boundary  in  the  ocean  is  important  during  ENSO  events. 

In  this  work,  the  forced  shallow  water  equations  are  solved  analytically  in  a  closed  equa- 


torial  ocr.  m  basin  using  a  method  of  solution  quite  similar  to  that  used  by  McCreary  (1976). 
Some  of  McCreary’s  simplifying  assumptions  were  dropped,  requiring  a  final  numerical  step 
to  construct  the  solutions,  and  a  wider  range  of  forcing  functions  was  investigated.  The 
problem  was  posed  as  an  initial  value  problem.  The  time  dependence  was  made  algebraic 
by  a  Laplace  transform.  The  dependent  variables  were  expanded  with  respect  to  latitude  in 
terms  of  the  Hermite  polynomials,  and  the  resulting  system  was  reduced  to  a  single  equa¬ 
tion  for  un(s,i),  the  transform  of  the  Hermite  coefficients  for  the  zonal  current.  Finally, 
the  solution  was  back-transformed  to  physical  space  numerically,  and  the  total  solution 
was  reconstructed. 

Section  2  reviews  the  system  of  equations  and  the  method  of  solution.  The  results 
obtained  using  two  different  forcing  functions  are  described  in  Section  3.  Some  final  com¬ 
ments  on  the  solutions  are  included  in  Section  4. 
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2  Basic  Formulation 


In  order  to  find  the  response  of  the  ocean  basin  to  prescribed  forcing,  the  shallow  water 
equations  on  the  equatorial  (3  plane  are  written  in  their  dimensional  form: 

dh  du 

9-^  + -  Pyv  =  T 

ox  dt 

dh  dv 

9-^+Pyu+~r~  =  0 

dy  dt 

dh  ,du  dv 

el  +  ddi+dai  =  0 

where  r  is  the  surface  wind  stress,  d  is  the  mean  depth  of  the  fluid,  and  the  other  variables 
have  their  usual  meanings.  The  independent  variables  are  nondimensionalized  according 


.  i  i  ~ 

/? 2  c 2  i 


/ 3*c  2x 


P*C  2y 


where  c  =  ( gd )  2.  The  fundamental  equations  for  the  system  are  then 


■hx  +  ut  —  yv  = 


a-1  -A 
p  2c  2r 


-hv  +  yu  +  vt 
c 


-ht  +  uz  +  Vy 
c 


The  zonal  current  is  required  to  vanish  at  the  ends  of  the  basin  (x  =  0,/),  and  the 


meridional  current  is  required  to  vanish  far  from  the  equator.  Initially,  the  system  is  at 
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rest,  with  a  flat  surface.  These  boundary  and  initial  conditions  are  expressed  as 


u(0,i,y) 

= 

0 

(4) 

u(0,x,y) 

= 

0 

(5) 

h(0,x,y) 

= 

0 

(6) 

u  —  0 

at 

x  —  0,1 

(7) 

v  -*  0 

as 

y  -»  ±<x> 

(8) 

The  values  of  the  constants  are  shown  in  Table  1. 


A  Laplace  transform  of  (l)-(3)  gives 


g  - 

-hx  +  su  —  yv 
c 

9 1  .  A 

-hy  +  yu  +  sv 


/T5<T*f 

0 


-sh,  +  ux  +  vu 
c 


=  0 


(9) 

(10) 

(11) 


where  |(s,x,y)  =  £(£)  =  /0°°  x,  y)e  3tdt.  h  is  eliminated  by  taking  s(9)  —  (ll)^  and 

*(10)  -  (H)v 


s2u  —  uxx  —  ysv  —  vXy  =  (3  *c  *sf  (12) 

ysu  -  uxy  +  s2v  -  Vyy  =  0  (13) 

The  dependent  variables  and  the  forcing  function  are  expressed  as  sums  of  the  Hermite 
polynomials,  Hn. 

oo 

i(s,x,y)  =  Cn(-s,x)Hn(y)e_y2/2 

n=0 
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dimensional  value  nondimensional  value 
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Expanding  (12)  and  (13)  in  terms  of  H„,  multiplying  by  Hm  e_y2/2  and  integrating  over  y, 
a  pair  of  equations  for  un  and  vn  are  obtained. 

(s2-D2)un-(n  +  l)(s  +  .D)vn+1-^(s  -D)vn^  =  (14) 

(n  +  l)(s  -  D) un+1  +  ~ (s  +  Z?)un_! 

—  (n  +  2)(n  +  l)5n+2  +  (>s2  +  n  +  ~)vn  -  =  0  (15) 

2  4 

These  equations  are  rewritten  in  matrix  form 

A|u)  —  B\v)  =  C\t) 

E\u)  +  Eju)  =  0 

where 

^n,m  —  [s  —  D  )6n,m 

Bn,m  =  ~(s  —  -D)£n,m+1  +  (n  +  l)(s  +  D)6nm-i 
Cn.m  =  P  *  C  *s6n  rn 

En,m  =  -($  +  F>)Snm+ 1  -h  (n  +  l)(s  —  D)6nm-i 

Fn.m  —  —  ~^n,m+2  +  (s2  +  n  -f  -)<fr»,m.  —  (n  +  2)(n  +  l)Sn  m-2 

and  D  -  The  matrix  operators  A,  B,  C,  E,  and  F  are  all  linear,  with  constant 
coefficients.  Then  FB  —  BF  =  0,  and 

{FA  +  BE)\u)  =  FC\t)  (16) 
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The  operators  in  (16)  are  expanded 


(. FA  +  BE )  =  -±(s2-D2)<5n,m+2  +s2  +  n  +  |(s2  —  D2)6nm  -(n  +  2)(n  +  l)(s2  -  D2)6n,m 

+  i (s*  -  D*)6nim+i  +\n(s  -  D)26n<m 

+  j  {ft  +  l)(-s  +  D)76nm  +(n  +  l)(n  +  2)(s2  —  D2)6nm 

FC  =  j3-*c~*s{-\6ntm+ 2  +  (s2  +  n  +  |)6„im  -(n  +  2)(n  +  l)6n,m_2} 

The  off  diagonal  terms  of  (FA  +  BE)  all  cancel,  so  that  (16)  may  be  written 

{s2Z?2  +  sD  +  (s2  +  2n  +  l)s2}  t2„  =  (3~*c~*s  j-^fn_2  +  (s2  +  n  +  ^)fn  -  (n  +  2)(n  +  l)fn+2 j 

(17) 

General  solutions  to  (17)  are  constructed  as  the  sum  of  two  parts,  un  =  Un  -f-  t)n, 

A 

where  Un  is  a  particular  solution  depending  on  the  forcing,  and  vn  is  a  solution  to  the 
homogeneous  form  of  the  equation.  Solutions  to  the  homogeneous  equation  are  given  by 
0n  =  fn(s)ekn(‘\  where  /„  is  an  arbitrary  function  of  s  (to  be  determined  by  the  boundary 
conditions),  and  kn(s)  is  given  by  the  dispersion  relation 

—s2k2  +  sk  +  (s2  -f  2n  +  l)s2  =  0 

k*  =  ~  -l±yjl  +  4s2(si  +  2n  +  l)  (18) 

Zo  L 

The  Hermite  representation  guarantees  that  the  solution  will  vanish  as  y  — *  ±oo.  The 
boundary  conditions  at  x  =  0,/  must  be  satisfied  by  an  appropriate  choice  of  fn(s)-  That 

a 

is,  the  fn(s )  must  be  chosen  such  that 

fn  (s)  +  fuis)  +  Un(s,0)  =  0 

ft(s)ek*1  +  f~(s)ek"1 +  Un(s,l)  =  0 
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Thus,  the  coeficients  /*(s)  are  given  by 


/*(*)  =  =F 


Un(sJ)-ekZ‘Un(s,  0) 


et+<  _ 


General  solutions  to  (17)  which  satisfy  both  initial  and  boundary  conditions  are  then 


found  for  different  forcing  functions  r(t,x,y),  and  back  transformed  to  physical  space. 


3  Solutions 


a  Zonally  uniform  forcing 


The  first  case  which  was  investigated  was  driven  by  a  simple  wind  stress  function 


r(t,  x,y)  =  —  (tO(t)  -  [t  -  t0)Q[t  -  t0)]  e  ^ 
to 


where  0  is  the  Heavyside  step  function.  (20)  represents  a  zonally  uniform  eastward  wind 
stress  with  a  Gaussian  meridional  profile.  The  width  of  the  Gaussian  corresponds  to  that  of 
a  first  meridional  mode  atmospheric  Kelvin  wave.  The  magnitude  of  the  forcing  increased 
from  zero  to  a  finite  value  corresponding  to  a  uniform  westerly  wind  anomaly  of  about  10 
m  s"1  over  a  period  of  about  one  month.  This  uniform  “spin-up”  forcing  is  very  similar  to 
that  used  by  McCreary  (1976). 

The  above  forcing  was  expanded  in  terms  of  Hermite  polynomials,  and  each  term  was 


Laplace  transformed  to  yield 


^  f1  -  «•-] T- 


where  T„  is  the  nth  Hermite  coefficient  of  At  *'"/2 .  Particular  solutions  to  (17)  are 


A  2  £  2  1  J  1 

Un(s’  X)  =  s3(s2  +  2n  +  1)  T0  [~4Tn-2  +  (s2  +  n+2)Tn~{n  +  2){n  +  l)Tn+2 


Usir.,,  (19),  the  general  solution  for  un  is  found  to  be 


z te - ^  (21) 


u 


The  Hermite  series  for  r(y)  is  truncated  at  n  =  6,  and  solutions  are  found  for  un(s,x). 
These  are  back  transformed  numerically,  using  the  IMSL  subroutine  INLAP,  and  the  so¬ 
lution  for  u(t,x,  y)  is  reconstructed  in  physical  space. 

The  zonal  current  anomaly  for  this  case  is  shown  in  Figs.  1  and  2.  These  figures 
show  an  initial  response  of  eastward  current  anomalies  along  the  equator  in  response  to 
the  eastward  stress  anomaly.  By  t  =  5  (roughly  8.5  days)  this  response  is  limited  to  a 
narrow  latitude  band  along  the  equator,  with  a  wider  region  of  eastward  flow  near  the 
eastern  end  of  the  basin.  The  strongest  response  at  this  time  is  the  .20  m  s'1  westward 
currents  located  on  the  equator  near  the  western  end  of  the  basin.  Less  intense  westward 
currents  flank  the  eastward  equatorial  flow  across  most  of  the  basin,  with  eastward  return 
flow  at  higher  latitudes. 

At  t  =  20  («  34  days)  the  forcing  has  reached  its  final  value.  The  equatorial  currents 
are  now  westward  at  all  longitudes,  with  a  maximum  of  .18  m  s_1across  most  of  the  basin. 
The  patch  of  enhanced  westward  flow  near  the  western  boundary  is  still  evident,  but  its 
magnitude  has  decreased  to  about  .14  ms-1.  Under  the  influence  of  continued  forcing,  the 
same  general  pattern  persists.  The  meridional  extent  of  the  westward  currents  continues  to 
expand  with  time  and  the  relative  magnitudes  of  the  eastern  and  western  current  maxima 
change  so  that  at  t  —  100  («  171  days)  they  are  both  «  1.6  m  s_1(Fig.  2). 

Only  the  far  western  part  of  the  basin  has  reached  a  steady  state  configuration  at  this 
time.  Examination  of  the  Hermite  coefficients  (not  shown)  indicates  that  only  the  n  =  0 
term  has  reached  a  steady  state  value.  The  n  =  2  term  has  equilibrated  over  the  western 
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U(X.Y),  T=20  (INT-0.02) 


. - 1  i - 1  : - 1  i  i  i  !  i 
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X 


Figure  1:  Zonal  current  anomaly,  a)  t  =  5,  contour  interval  =  .002  m  s_1;  b)  t  =  20, 

contour  interval  =  .02  ms  1 .  (Westward  currents  are  dashed.) 
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Figure  2:  Zonal  current  anomaly,  a)  t  =  50,  contour  interval  =  .1  m  s-1;  b)  t  =  100, 


contour  interval  =  .2  m  s_1. 
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one  fourth  of  the  basin,  and  higher  order  terms  over  successively  smaller  regions  near  the 


western  boundary.  This  is  to  be  expected,  because  of  the  slower  phase  speeds  associated 
with  the  latter  waves.  However,  since  the  magnitudes  of  the  lowest  order  coefficients 
are  much  greater  than  those  of  the  higher  order  terms,  the  qualitative  character  of  the 
solution  does  not  change  significantly.  The  steady,  basin-wide  eastward  wind  stress  results 
in  a  relative  rise  in  sea  level  in  the  eastern  part  of  the  basin,  and  a  westward  equatorial 
flow  with  weak  eastward  return  flow  at  higher  latitudes. 


b  Eastward  moving  pulse 

The  second  case  which  was  studied  is  more  realistic,  and  therefore  somewhat  more  complex. 
The  wind  stress  is  given  by 


r(t,  x,y)  —  [l  -  e  at  [0(x  -  gt)  —  0(x  -  A  —  ryt)]  At ~7^7 


Thus,  the  surface  stress  is  assumed  to  be  in  the  form  of  an  eastward  moving  pulse  (top-hat) 
of  anomalous  westerly  winds  similar  to  that  used  by  McCreary  and  Lukas  (1936)  in  their 
study  of  an  unbounded  equatorial  ocean.  The  zonal  extent  of  the  pulse  (A)  is  taken  to 
be  half  the  length  of  the  basin  («  7000  km)  and  it  moves  eastward  with  speed  g  such 
that  it  traverses  the  basin  in  about  16  months.  Its  magnitude  grows  from  zero  to  its 
asymptotic  value  with  an  exponential  time  constant  of  about  four  months.  During  this 
time,  the  leading  edge  of  the  pulse  travels  from  the  western  boundary  of  the  ocean  to 
x  ss  12.  The  meridional  shape  and  final  magnitude  of  the  forcing  are  the  same  as  in  the 
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zonally  uniform  case  discussed  above.  This  wind  stress  anomaly  is  more  realistic  than  the 
uniform  stress  discussed  in  the  previous  section  in  that  it  is  not  zonally  uniform,  and  the 
region  of  strongest  anomalous  wind  stress  propagates  towards  the  east,  eventually  leaving 
the  region.  This  behavior  is  not  unlike  the  wind  anomalies  that  are  observed  during  major 
ENSO  events  such  as  1982-83,  and  it  is  similar  to  that  used  by  Blundell  and  Gill  (1983) 
and  McCreary  and  Lukas  (1986). 

The  Laplace  transform  of  the  nth  term  of  the  Hermite  expansion  for  r(t,x,y)  is 

fn(s,x)  =  Q{l  —  e-Wn_0(x_A)[i_e-.(*-A)/,J]} 

{l  -  e-(4+a)l/r’  -  0(x  -  A)  [l  -  Tn 


A  particular  solution  to  (17)  with  this  forcing  is 

Un(s,x)  =  An(s){pn(s)  +  Q„(s)e-,I/'7  +  JR„(s)e-(,+  °)I/r' 

-0(x  -  A)  [Pn{s)  +  Qn{s)e-’{z~^/ri  +  i2n(s)e-^+a)(l-A^'7]  } 

where 


A„(s) 

Pn(s) 

Qn(s) 

Rn{s) 


0  7C  7  —  4^-2  +  (-s2  +  n  +  -)T„  -  (n  +  2)(n  +  l) T^+2 j 
a 

s2(s  +  a)(s 2  +  2n  +  1) 

_ _ 

s2[r)2(s2  +  2n  +  1)  -  tj  -  s2] 

_ rf _ 

(s  +  a)[sr}2(s2  +  2n  +  1)  —  (s  +  a)y  —  s(s  +  a)2] 
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Since  Un(s,x)  is  discontinuous  at  x  —  A,  the  functions  /*($),  discussed  in  Section  2,  must 
be  modified  slightly  to  allow  different  solutions  for  x  <  A  and  x  >  A.  Additional  boundary 
conditions  require  that  the  total  solution,  u„(.s,x),  and  its  first  derivative  with  respect  to 
x  are  continuous  at  x  =  A. 

The  Hermite  coefficients  of  the  solution,  u„(t,x),  are  shown  in  Figs.  (3)-(7).  Note 
that,  in  order  to  show  details  of  the  response,  each  figure  has  a  different  scale  on  the 
ordinate.  The  points  at  which  the  curves  appear  to  be  jagged  (e.  g.  u0(50,x)  for  x  ss  20) 
correspond  to  points  at  which  the  numerical  inversion  of  the  Laplace  transform  failed  to 
converge  to  the  desired  (1%)  tolerance,  and  a  10%  tolerance  was  accepted. 

The  dominance  of  the  lowest  order  (Kelvin  wave)  term  is  clearly  seen.  Immediately  after 
the  pulse  of  eastward  wind  stress  is  initiated  the  n  =  0  response  retains  the  shape  of  the 
forcing,  but  propagates  eastward  at  a  considerably  greater  speed  than  does  the  prescribed 
forcing.  This  is  the  result  of  a  superposition  of  a  series  of  Kelvin  waves,  each  propagating 
with  speed  c,  which  is  about  5.6  times  greater  than  the  speed  of  the  forcing,  r?.  Once  the 
leading  edge  of  the  Kelvin  wave  front  reaches  the  eastern  boundary  (t  «  25  «  43  days)  the 
shape  of  the  pulse  is  no  longer  preserved,  due  to  reflections.  The  strongest  response  is  seen 
at  t  =  100  (zz  170  days).  By  t  =  200,  the  leading  edge  of  the  forcing  has  moved  beyond  the 
eastern  boundary  of  the  ocean  basin,  so  that  only  the  eastern  third  of  the  basin  is  subject 
to  direct  forcing. 

The  n  =  2  (Fig.  4)  response  is  similar  to  that  for  n  =  0,  except  that  the  leading  edge 
of  the  response  is  retarded,  due  to  the  slower  phase  speed  of  the  higher  order  wave  mode. 
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The  magnitude  of  the  response  is  much  smaller  because  the  direct  forcing  of  this  mode  is 
also  much  smaller.  The  higher  order  response  (Figs.  5-7)  follow  this  same  pattern. 

Even  though  the  magnitude  of  the  n  =  0  term  is  considerably  greater  than  that  of  the 
higher  order  modes,  the  latter  have  a  distinct  influence  on  the  character  of  the  complete 
solution.  Fig.  8  shows  the  total  zonal  current  anomaly  at  t  =  5  and  t  —  50.  As  might  be 
expected  the  initial  response  of  the  western  part  of  the  basin  is  rather  similar  to  that  found 
for  the  uniform  forcing  case.  The  currents  are  weaker,  because  the  forcing  approaches  its 
final  value  more  slowly,  but  the  basic  structure  of  westward  currents  at  at  the  equator, 
flanked  by  eastward  return  flow,  is  the  same.  At  t  =  5,  the  leading  edge  of  the  forcing 
is  at  x  =  25.9.  The  Kelvin  wave  front  is  at  x  =  30.  To  the  east  of  i  =  30  there  is  no 
disturbance,  as  the  signal  has  not  yet  reached  this  region.  By  t  =  50,  the  leading  edge  of 
the  forcing  has  reached  x  «  35.  To  the  west  of  this  point,  the  solution  still  resembles  that 
of  the  earlier  case.  However,  evidence  of  the  trailing  edge  of  the  forcing  begins  to  appear 
in  the  far  western  part  of  the  basin.  Just  ahead  of  the  leading  edge  of  the  wind,  a  burst  of 
eastward  motion  (a  Kelvin  wave  front)  exists,  with  only  a  weak  response  near  the  eastern 
boundary. 

At  later  times  (Fig.  9)  the  effect  of  the  nonuniform,  transient  nature  of  the  forcing 
becomes  more  pronounced.  At  t  —  100  moderate  («lm  s-1)  westward  currents  dominate 
the  center  of  the  basin.  These  currents  decrease  more  rapidly  towards  the  west  than  in  the 
previous  case  (see  Fig.  2)  because  the  forcing  there  has  vanised.  The  Kelvin  wave  front  of 
strong  eastward  currents  ahead  of  the  forcing  pulse  has  gained  strength,  and  now  extends 
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a  considerable  distance  in  the  meridional  direction.  By  t  =  200  (ss  11  months),  only  the 
eastern  third  of  the  basin  is  directly  influenced  by  the  forcing.  There,  the  strong  westward 
currents  persist  near  the  equator,  with  eastward  flow  to  the  north  and  south.  To  the  west 
of  the  active  forcing  region,  the  flow  is  weaker. 

At  still  later  times  (Fig.  10)  the  pulse  of  anomalous  wind  stress  has  exited  the  model 
ocean  basin.  For  t  >  380,  the  forcing  vanishes  across  the  entire  basin.  The  current 
anomalies  which  persist  are  due  to  reflections  of  the  waves  which  were  generated  as  the 
forcing  traversed  the  basin.  Because  the  model  includes  no  dissipation,  these  waves  will 
persist  for  all  time,  but  the  shape  of  the  disturbance  will  continue  to  fluctuate  as  the 
different  wave  modes  superpose  in  different  ways.  At  t  =  400  the  region  of  westward  flow 
has  expanded  to  fill  almost  the  entire  basin;  the  familiar  bands  of  eastward  flow  have  all 
but  vanished.  While  this  is  not  a  strictly  steady  state  configuration,  the  genera!  state 
persists  through  t  =  500,  with  only  small  qualitative  changes  due  to  the  slow  phase  speeds 
of  the  higher  order  modes  which  dominate  at  this  time. 
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4  Sumary  and  Conclusions 


In  this  report  we  have  presented  the  analytic  solution  to  the  shallow  water  equations  in  a 
closed  basin  subject  to  arbitrary  wind  forcing.  Two  physically  realistic  example  solutions 
were  shown  and  interpreted  in  some  detail.  One  example  corresponded  to  the  spin  up  of 
an  equatorial  ocean  basin  to  a  zonally  uniform  westerly  wind  anomaly.  The  other  example 
was  the  transient  response  of  the  basin  to  a  zonally  oriented  top-hat  wind  pulse  that  moved 
across  the  basin  in  a  period  of  16  months.  This  represented  the  wind  forcing  of  an  idealized 
version  of  the  atmospheric  Southern  Oscillation. 

The  analytic  solution  resembles  published  results  for  an  ocean  basin  (e.  g.  Cane  (1979); 
Busalacchi  and  O’Brien  (1981);  Weisburg  and  Tang  (1983))  which  until  now  were  obtained 
only  using  numerical  models.  Although  our  results  are  somewhat  incomplete  without  the 
rest  of  the  solution  (i.  e.  v  and  h),  it  is  clear  that  the  present  analytic  solution  offers  con¬ 
siderable  potential  for  interpreting  observations  and  numerical  model  results  in  equatorial 
oceanography. 
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